The Universal Gaussian in Soliton Tails 
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We show that in a large class of equations, solitons formed from generic initial conditions do not have 
infinitely long exponential tails, but are truncated by a region of Gaussian decay. This phenomenon 
makes it possible to treat solitons as localized, individual objects. For the case of the KdV equation, 
we show how the Gaussian decay emerges in the inverse scattering formalism. 



Recently, progress has been made understanding the 
time development of the leading exponential edge in 
propagating fronts JnJ . Inspired by results on fronts in the 
Fisher equation J2|,||], it was shown in ffl that in generic 
reaction-diffusion equations like Ginzburg-Landau (GL), 
an initial condition with a compact front gives rise to 
a front with a leading exponential edge that does not 
extend forever, but rather is cut off a finite distance 
ahead of the front. The need for this is clear on physical 
grounds: the front, if one is sufficiently far ahead of it, has 
not had time to make its presence felt, and so the field ex- 
hibits the typical Gaussian falloff of the Green's function 
(which lacks an intrinsic scale). It was discovered that 
there is a well-defined transition region of width O(Vi) 
wherein the field crosses over from the steady-state ex- 
ponential to Gaussian falloff. This "precursor" transition 
region propagates out ahead of the front, with a velocity 
c* which is greater than twice the velocity c of the front 
itself. 

Given the basic underlying physics, the existence of 
such transition regions must be a very general phe- 
nomenon, true not just of reaction-diffusion fronts, but 
of other propagating solutions, such as the solitons in the 
KdV equation, 

Ut = ~ u xxx + Guu x . (1) 

The known exact one-soliton solutions take the form 
u(x,t) = —\c sech 2 ( | ^fc(x — xo — ct) ) , and exhibit expo- 
nential decay for large x. The inverse scattering trans- 
form 0] tells us that the solution of KdV with generic 
initial condition u(x, 0) = 4>(x) (with <fi(x) — > suffi- 
ciently rapidly as |a;| — > oo) consists of a train of soli- 
tons moving to the right, along with a dispersive wave 
travelling to the left. As above, we can argue that when 
4>{x) has compact support, the solution emerges from the 
rightmost soliton decaying as exp(— \fcx), where c is the 
relevant speed, but for sufficiently large x the presence 
of the solitons will not yet be felt, and the behavior of 
the solution will be determined by the Green's function 
of the linearized equation ut — —u xxx , i.e. u will decay 
roughly as exp(— 2ir 3 / 2 /3v / 3t) ||. Thus there is a tran- 
sition in the nature of the decay. How and where does 



this transition take place? 

If there is more than one soliton in the soliton train, 
say two, with speeds c±, c 2 (c 2 > Ci > 0), then a further 
problem arises. The solution emerges from the faster- 
moving soliton decaying as exp(— ^foix), but because the 
tail of the slower-moving soliton falls slower, it is possible 
that it will return to dominate, i.e. the decay will slow 
to exp(— yfcix). (For very large explained above, 

the solution must go roughly as exp(— 2x 3 / 2 /3\/3t)-) We 
note that exact two-soliton solutions 

. . (c 2 - ci)(c 2 cosh 2 a\ + ci sinh 2 a 2 ) 
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where a\ = h-y/c~i(x — c\t) and a 2 = \\fc^{x — c 2 i), 
exhibit exactly this phenomenon: when x is large u ~ 
exp(— y/cix), i.e. the tail of the slow soliton dominates 
the decay. This is clearly undesirable in physical situa- 
tions, as it means solitons cannot be considered as iso- 
lated objects. 

In this paper we study the tails of generic solitons, i.e. 
those produced by taking a generic compactly supported 
(or very rapidly decaying) initial condition for KdV. The 
results are quite striking. If the soliton is centered at 
x = ct, then in a region of width O^IH 1 / 2 ) around 
x = c*t, where c* = 3c, there is a rapid transition in 
the behavior of the tail, and the decay changes from 
exponential to Gaussian. This behavior persists until 
x — c*t = 0{c 1 l 2 t 2 / 3 '), and then there is a second tran- 
sition to the final region, in which u decays roughly as 
exp(— 2cc 3/ ' 2 /3v / 3£)- In the Gaussian region the decay is 
very rapid, and the soliton tail is effectively cut off, mak- 
ing the soliton an isolated object. 

This behavior is predicted by analysis of the linearized 
equation u t — —u xxx alone, and does not involve any 
of the special properties of the KdV equation. It is thus 
universal for soliton solutions of PDEs with this lineariza- 
tion, and the existence of a Gaussian cut-off region is in 
fact universal in a much larger class of equations, and is 
responsible, as we have explained, for the individuality 
of solitons. The existence of exact two-soliton solutions 
in which the solitons are not genuinely separate objects 
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is one of the special properties of the KdV equation, and 
is not physical. 

The rest of this paper proceeds as follows: We first give 
the arguments for the behavior of the tail outlined above. 
We then show how this makes the soliton an isolated 
object. In the last part, we connect our results to the 
exact Inverse Scattering solution, and use this to generate 
a numerical solution of the KdV equation, confirming the 
picture presented. 

Soliton Tails — Since we only want to look at the soli- 
ton tail, where u is small, we work with the linearized 
equation u t = —u xxx . We have explained in the introduc- 
tion why the standard soliton tail u — — 2cexp(— y/c(x — 
ct)) is not a physically acceptable solution to this; we 
require a solution that for large x decays faster. We look 
for an acceptable solution in the form 

u(x, t) = ~2cexp(~^/c(x — ct))f(y, t) (3) 

where y = x — c*t and c* is a constant to be determined, 
and where / satisfies boundary conditions / — > 1 as y — > 
— oo and / — > as y — > +oo. Substituting in, we find / 
satisfies 

ft = -fyyy + ?>J~ c fyy + ip* - 3c)f y . (4) 

We choose c* = 3c so that the term in f y drops out. 
The key point is that at large times the diffusive term, 
3y/cfy y , dominates the RHS out to y of order c 1 / 2 ^/ 3 , 
whereas the dispersive term, f yyy , dominates for y 3> ct. 
To see this, let us first drop the f yyy term. The resulting 
diffusion equation then has the exact scaling solution 

f(y, t) = -erfc ( —jJL -) (5) 

yy ' 2 V2V3c 1 /^i/2; v ; 

which can be seen to satisfy the boundary conditions. 
This solution transforms the original exponential falloff 
of u to a Gaussian falloff in a region of width Oi^'H 1 / 2 ) 
around y — 0, or equivalently x = 3ct. This erfc cutoff of 
the exponential moving out ahead of the front is exactly 
what occurs in the case of the GL equation there it 
is an exact solution of the relevant linearized equation. 

Here, however, we have to address the effect of the ne- 
glected f yyy term. There are two cases, depending on the 
size of the scaling variable z = y/c 1 / 4 ^ 1 / 2 . For z — O(l), 
the f yyy term is down by a factor c -3 / 4 ^ 1 / 2 and so is in 
fact negligible for large t ^> c~ 3 l 2 . For large z, on the 
other hand, the f yyy term induces a correction of order 
z 3 jc 3 ^t 1 ! 2 and so once z = 0(c 1 / 4 t 1 / 6 ), or equivalently 
y = 0(c 1 / 2 t 2 / 3 ), it is no longer negligible, as we noted 
above. In fact, for very large y 3> ct, the f yyy dominates, 
and the controlling factor in / is the exp(— y 3 / 2 /t 1 / 2 ) of 
the Green's function, as expected. In this regime, the 
diffusive term can be seen to be of sublcading order, in- 
ducing a correction to the argument of the exponential 
of order \J~cy. 



The upshot of this is that the cutoff of the exponential 
is provided by the er/c, inducing Gaussian decay. Only 
much later does the Gaussian decay slow down to that 
prescribed by the Green's function. It is clear from the 
very general nature of this argument that this scenario is 
universal, applying to a wide range of soliton equations, 
as well as front propagation problems like GL. 

We have shown that the erfc solution is consistent, but 
because the basic equation (||) is linear we can go further 
and prove that it in fact is what arises from solving the 
initial value problem for compact initial conditions. To 
do this, we solve (Q) by taking a Fourier transform in 
space, finding 

f(y,t) = — e ik «-^ ck2 - ik3 »f(k)dk, (6) 

^ n J—ca+ie 

where f(k) is the Fourier transform of the initial condi- 
tion: 

/oo 
f{y,Q)e- ik Hy . (7) 
-oo 

It will turn out that in the region of y <C ct, the only 
thing we need to know about / is its behavior at small 
k. Since the long-distance structure of f(y, 0) is that of a 
step-function 9(—y), the small k behavior of / is precisely 
that of the Fourier transform of 0(—y), i.e. for small |fc| 

/>) w \ , Im(fe) > . (8) 

The range of validity of this approximation is |fc| <C a -1 , 
where a is a characteristic length scale of f(y, 0). 

We now use saddle point techniques to evaluate the in- 
tegral ^ . We denote the factor in the exponential in the 
integrand by g{k), i.e. g{k) = iky — (3^/ck 2 — ik 3 )t. For 
all positive t and y > —Set, g{k) has two critical points 
on the imaginary axis, at k± = i(±\/c + y/3t — y/c). For 
y > we deform the integral in (0) to the steepest de- 
scent integral coming in along the ray arg(fc) = 57r/6, 
going through the critical point on the positive imag- 
inary axis at k = k + , and going out along the ray 
arg(fe) — tt/6 (along both these rays ik 3 is negative real). 
Writing A = (3xt) 1 ^ 4 (k — k + ), the factor in the expo- 
nential is g(k) = g(k + ) - A 2 + i\ 3 (t/27x 3 ) 1/4 ; provided 
t/x 3 is small, which it will be, for example, for posi- 
tive y and t 3> c~ 3 / 2 , we can ignore the term in A 3 . 
Turning now to the factor multiplying the exponential 
in @, this is f{k) = f{k+ + (3xt)- 1 / 4 A). Provided |fc + | 
and (3xi) -1 / 4 are sufficiently small (y <C min(ct, y/ca^ x t) 
and t a 2 c~ 1 / 2 are sufficient conditions), we can use 
the approximation (|[) to estimate this. Putting all 
the approximations together, we have that for y > 0, 
y <C min(c£, y / ca~ 1 t) and t > max(a 2 c -1 / 2 , c~ 3 / 2 ): 

1 r°° ie~ x2 
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[g(k + ) + (3xt)^ 2 \k + \'- 
erfc((3:ci) 1/4 |fcH 



2 6XP 



(9) 



(See for example M 3.466 for the necessary integral.) 
This expression can be substantially simplified provided 
y <C 6c 1 / 2 i 2 / 3 , in which case the factor in the exponential 
m (g) is o(l), and we recover (ph. A similar calculation 
recovers (||) in the case y < as well. 

We can also use the Fourier integral @ to learn more 
about the post-Gaussian region. As y gets larger, so does 
|fc+| and we can no longer use the small- k limit of /. The 
physical meaning of this is that this region is sensitive to 
the exact initial conditions. This must be so, since the 
large y regions feel only the initial condition, since no 
other information has had time to propagate there. Nev- 
ertheless, we can obtain the controlling factors in f(y, t). 
Once we are sufficiently far from the k — pole in /, we 
can treat / as a constant when doing the saddle point 
integral. Using the exact form of k + above we obtain 



f(y,t) 



27T 1 /2(3te)l/' 



■ exp 



2x fx 

T V zt 



+ ^fc{x - ct^j 
(10) 



This equation is valid for in the regime t 3> max(a 2 c 1 I 2 , 
c~ 3 / 2 ) and z» 1, and also in the regime y ^> max(a 4 t _1 , 
c 1/4 i 1/2 ), t arbitrary. 

We see that the two approximations, Q and ( |l0| ) have 
a region of overlap for large t, namely cr^i 1 / 2 « j « 
c 1 / 2 t 2 / 3 . This is precisely the Gaussian region. Either by 
using y <C ct and the pole approximation for / in ([l(]) or 
by using the large argument approximation of erfc in (j^) 
we recover 



f(y,t) 



• exp 



y 



12^/ct 



(ii) 



Soliton Individuality — As explained in the introduc- 
tion, if solitons had infinitely long exponential tails, then 
in a situation where there were two solitons, with speeds 
ci,C2 (c2 > ci > 0), the tail of the slow soliton would 
dominate the large x decay. In this section we show that 
the cutting off the tail of a soliton with speed c at x = 3ct 
(up to an additive constant) prevents this. 

Assuming exponential tails, the point x(t) where the 
slow soliton tail returns to dominate is determined by an 
equation of the form ^/c\(x{t) — c\t) — ^fo2.(x{t) — c^t) + 
cnst and so 



The meaning of this, as emphasized in the introduc- 
tion, is that solitons can genuinely be regarded as isolated 
objects; this is only due to the Gaussian cut-off behavior 
we have discussed. 

Connection to Inverse Scattering — The KdV equation 
can be solved via the Inverse Scatting Transform (1ST). 
This should reproduce the results obtained by our general 
arguments above. For large x, where the field u(x,t) is 
small, the 1ST tells us that u rj 2-^F(x,t), where 



F(*,t) = 5> 2 e 8 ^- 2 -* + ij Kt) 



e 8ik 3 t+2ikx ( y c 



(13) 



The various constants in this equation are scattering data 
for the Schrodinger operator — d 2 4- <fr(x) associated with 
the initial condition 4>{x)\ specifically — k 2 gives the dis- 
crete spectrum (n = 1, . . . , N), the c n are the normaliza- 
tion constants and b(k) is the reflection coefficient. 

When 4>(x) is a reflectionless potential, b(k) vanishes 
and exact soliton solutions are obtained. However, for 
generic initial conditions, including the compact initial 
conditions considered here, b(k) does not vanish. In 
fact, it has a number of poles on the positive imaginary 
axis, which are associated with the solitons that develop. 
When we move the integration path in (|l^) to a new con- 
tour C above all these poles, the residues exactly cancel 
the discrete sum in (|l3]) |1, leaving, after a rescaling of 
k by a factor 2, 

u(x,t) = — I exp(ln6(fc/2) + ik 3 t + ikx)kdk , (14) 
2tt Jc 

Translating k by 2iki, where in\ is the location of the up- 
permost pole of b(k) (which gives rise to the fastest soli- 
ton, with velocity c — 4k 2 ), we obtain the solution in ex- 
actly the form given by equations (||) and (||) , where / (k) 
is identified with (Ki~ik/2)b(k/2 + iKi)/c. As expected, 
this f(k) has a pole at zero, and no other singularities 
in the upper-half fc-plane. The only trivial difference is 
that the residue at the pole at zero is not i, but some 
multiple thereof, which lets us calculate the phase shift 
of the soliton, information not accessible in the general 
framework. 

We can actually use ([l4]) to numerically calculate u 
for a given initial condition, and thus verify our analytic 
arguments. We present results for the case where <f)(x) is 
a square well of depth V and width 2a centered around 
x = 0: for this case we have 



dx(t) 
dt 



3/2 _ 3/2 



Ve 



-2iak 



1/2 



J/ 2 



= ci + v^+c 2 . (12) m = {q+ikcotqa)(q _ iktimqa) • i = vvn 



Thus the putative point of "return of the slow soliton" 
must move with a speed c\ + ^/ c\Ci~\- ci. But this exceeds 
c\ = 3ci , so for sufficiently large time the slow soliton has 
been cut off before it can return to dominate the decay. 



(15) 

The number of bound states (solitons) is 1 + [2^/Va/Tr], 
associated with poles on the imaginary axis of b(k). The 
poles all have |fc| < \fV . 



3 



1.0 



0.8 



£N 0.6 



z 04 

f 0.2 



0.0 



- 


»0»»MP00> 


o t=10 






o t=100 


- 




□ t=400 






1 ^OOJOO i. L=-S-S-?5jgW?F 



-150 -100 



-50 





x-3ct 



50 



100 



150 



FIG. 1. f(y,t) = -u{x,t) exp(c 1/2 (x - K - rf))/2c vs. 
y — x — Set for t = 10 (circle), 100 (diamond) and 
400 (square), starting from the square-well initial condition, 
u(x,0) — —V8(a — \x\), with V = o = 1. The solid line is the 
analytic approximation, (^). 

Separating the exponent in ( |l4| ) into its real and imagi- 
nary parts, A(k)+i&(k), the idea is to integrate along the 
contour of constant-phase $ = passing through the sad- 
dle point which lies on the imaginary axis. The integrand 
is symmetric about the imaginary axis, and is strictly 
decreasing as we move away from the saddle. Denot- 
ing our steepest-descent contour by fc(s) = t(s) + iuj(s), 
parametrized by the arclength s, u is given by u(x, t) — 
4X(oo)/7r, where 



Z(s) = f ds' L^- + r^j expos')) 



(16) 



We simultaneously find this steepest-descent contour and 
the integral I(s) along the contour by solving via Runge- 
Kutta the following third-order system: 



T = -Qjy/^ + Ql- A0>$ r /(a>2 + $2 ■ 

Z = + $ r r) exp(^)/v/^ + $2_ 



(17) 



Here the dot denotes a derivative with respect to s, and 
A is an arbitrary (positive) Lagrange multiplier parame- 
ter stabilizing the <I> = constraint. It is easy to verify 
that $(fc(s)) = — A$. In practice, Z decays rapidly in 
s, and the integration may be halted after sufficient ac- 
curacy is achieved. This process is easily repeated for 
various x, t, yielding the results in Figs. 1 and 2. In Fig. 
1, we plot f(y,t) = — u(x, t)/2cexp(— \fc{x — xq — ct)) 
for t = 10,100,400. We also exhibit our universal ap- 
proximation for /, (|B|). The argument is good for the 
smaller t, and excellent for the larger t's. In Fig. 2, we 
again plot f(y,t), this time in semi-log scale, along with 
our approximations (Js[) and (|l(]). We see that the er/c 
works some way past the exponential cutoff, and that 



( |Io| ) works from the Gaussian regime outward, in accord 
with our analytical understanding. 

In summary, we have found that in the KdV equation, 
and many other systems possessing solitons, the leading 
exponential edge of the soliton is only built up over time. 
At any given time, it is cut off at some point by a mul- 
tiplicative erfc factor, which transforms the exponential 
decay into a Gaussian falloff. The cutoff point moves 
out with time, so that the length of exponential edge 
increases linearly with time. This cutoff phenomenon 
serves to give the soliton an individual identity. 
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FIG. 2. Same data as in Fig. 1, in semi-log scale,plotted 
together with the two overlapping approximations, (H) (solid 
line) and (Q (dotted line). 
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